Entangling two distant nanocavities via a waveguide 
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In this paper, we investigate the generation of continuous variable entanglement between two 
spatially-separate nanocavities mediated by a coupled resonator optical waveguide in photonic crys- 
tals. By solving the exact dynamics of the cavity system coupled to the waveguide, the entanglement 
and purity of the two-mode cavity state are discussed in detail for the initially separated squeezing 
inputs. It is found that the stable and pure entangled state of the two distant nanocavities can 
be achieved with the requirement of only a weak cavity-waveguide coupling when the cavities are 
resonant with the band center of the waveguide. The strong couplings between the cavities and the 
waveguide lead to the entanglement sudden death and sudden birth. When the frequencies of the 
cavities lie outside the band of the waveguide, the waveguide-induced cross frequency shift between 
the cavities can optimize the achievable entanglement. It is also shown that the entanglement can 
be easily manipulated through the changes of the cavity frequencies within the waveguide band. 

PACS numbers: 42.50.Dv; 03.65.Yz 
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I. INTRODUCTION 

Entanglement of electromagnetic field with continuous 
variables [1], which embodies quantum correlations in 
amplitude and phase quadratures of the field, has been 
proven to be very important for building up continu- 
ous variable quantum communication network and 
quantum computation ^-^8^ . Generally, continuous vari- 
able entanglement is generated via nondegenerate optical 
parametric oscillation 0, [^, [9l-[T]| or beam splitters incor- 
porating with squeezed light beams 2, 12]. 

Based on recent achievements in the light manipula- 
tion with highly controllable photonic crystals we 
investigate in this paper the generation of continuous 
variable entanglement between two spatially-separated 
nanocavities mediated by a coupled-resonator optical 
waveguide (CROW) in photonic crystals. As we know, 
the long-distance entanglement is an essential ingredient 
for transmitting quantum information over long distances 
in quantum communication networks [3, EBl • Very re- 
cently, the distant entanglement between two qubits me- 
diated by a plasmic waveguide or two vibrating 
trapped ions assisted by a phononic reservoir (l7j has 
just been studied. Here we shall focus on the distant 
continuous variable entanglement between two nanocav- 
ity fields. The main merits of the present scheme are as 
follows. The nanocavity could be a point defect created 
in photonic crystals and its frequency can be simply con- 
trolled just by changing the size or the shape of the defect 
p^ . While the waveguide, as the mediator of the two 
cavities as well as the output channel, can be considered 
as a set of linearly coupled defects in photonic crystals in 
which light propagates due to the coupling between the 
adjacent defects [3,123]. The transmission properties of 
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the CROW can also be easily manipulated by changing 
the modes of the resonators and the coupling configura- 
tion [2l|. Furthermore, the coupling of the cavity to the 
waveguide is also controllable through the change of the 
distance between the corresponding defects 22] . Besides 
the above flexible controllability, the possibility of minia- 
turizing the entanglement setup with solid-state photonic 
structures is highly desirable for scalable and on-chip 
photonic quantum information processing [23l . [2^ . Obvi- 
ously, the bulk setups for producing continuous variable 
entangled light in Refs.jl, H, [5l-[T2j are not suitable for 
integration. In addition, compared to qubit or phononic 
entanglement proposed in Refs.p^ [T7|. optical entangle- 
ment with continuous variables can be easily detected 
and manipulated experimentally, which can make it very 
promising for the implementation of continuous variable 
quantum protocols 

The effects of various environments on continuous vari- 
able entanglement of optical fields or harmonic oscillators 
have been extensively investigated, where one is mainly 
interested in the non-Markovian decoherence dynamics 
of the entanglement [25- 29,]. For the present system 
of two spatially-separated nanocavities mediated by a 
waveguide, the waveguide actually serves as a structured 
reservoir which is highly controllable in experiments. We 
therefore are interested in the entanglement generation 
and its manipulation through the controls of the waveg- 
uide as well as the nanocavity properties. We shall study 
the exact entanglement dynamics solved from the exact 
master equation of the nanocavity system coupled to the 
waveguides that we developed very recently |30l433| . The 
temporal behavior of the entanglement and purity of the 
two-mode cavity state is discussed in detail for initially 
separated squeezing inputs. We find that when the cav- 
ity frequencies are resonant with the band center of the 
waveguide, a stable and pure entangled state of the two 
distant cavities can be generated with a requirement of 
only a very weak coupling between the cavities and the 
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FIG. 1: A schematic plot of two spatially-separated nanocav- 
ities coupled to a CROW structure at the sites rii {i — 1,2) 
in photonic crystals. 



waveguide. It shows that the strong cavity- waveguide 
coupling will lead to the entanglement sudden death and 
sudden birth (sj. When the frequencies of the cavities 
locate outside the waveguide band, the cross frequency 
shift of the nanocavities induced by the waveguide can 
optimize the achievable entanglement. In addition, it 
is also shown that the entanglement can be easily ma- 
nipulated by changing the cavity frequencies within the 
waveguide band. 

The rest of the paper is organized as follows. In Sec. II, 
we formulate the entanglement dynamics of two spatially- 
separated nanocavities coupled to a waveguide within the 
framework of the exact master equation and the correla- 
tion matrix, the later determines the continuous variable 
entanglement of the two nanocavities through the mea- 
sure of logarithmic negativity. In Sec. Ill, the properties 
of the entanglement and the corresponding purity are in- 
vestigated in detail. Finally, a conclusion is given in Sec. 
IV. 



II. HAMILTONIAN AND ENTANGLEMENT 
MEASURE 

As schematically shown in Fig. [1] we consider two 
single-mode nanocavities coupled with a coupled res- 
onator optical waveguide (CROW) in photonic crystals 
at different sites. Each cavity is formed by a point defect 
created in photonic crystals, with the cavity frequency 
tunable by changing the geometrical parameters of the 
defect. The waveguide in photonic crystals consists of a 
series of coupled point defects in which light propagates 
due to the coupling between the adjacent defects. Ex- 
perimentally, the large-scale CROW consisting of more 
than one hundred coupled resonators has been success- 
fully fabricated jlSj]. To be specific, let the cavity 1 be 
coupled to the waveguide at the site ni and the cav- 
ity 2 coupled to the waveguide at the site n2, as shown 
in Fig. [TJ By treating the waveguide as a tight-binding 
model, the Hamiltonian of the whole system is given by 

H 



H 



1=1,2 



k i,k 



where 

Wfc = Wo - 2^0 cos k, V^k = ^ sm{n,k), (2) 

with < fc < TT. The operators and are the anni- 
hilation and creation operators of the cavity fields with 
frequencies w^. The annihilation and creation operators 
bk and 6^, describe the Bloch modes uJk of the waveguide 
with ujQ being the identical frequency of each resonator in 
the waveguide. The frequencies tUi and ojq are tunable by 
adjusting the geometrical parameters of the correspond- 
ing defects. The strength characterizes the photon 
hopping between two adjacent resonators in the waveg- 
uide and is controllable by changing the corresponding 
distance between the two defects. The controllable cou- 
pling strength is the coupling of the ith cavity to the 
waveguide at the sites n^. We should point out that the 
frequencies of the two cavities and the waveguide band 
considered in the above system should lie inside the pho- 
tonic band gap of the photonic crystals. Then the photon 
loss into the photonic crystals is totally negligible [35|. 

We will investigate the generation of the entanglement 
between the two spatially-separated cavity fields through 
the controllable waveguide, where the two nanocavities 
are initially prepared in a separated two-mode Gaussian 
state. For a two-mode Gaussian state, its quantum statis- 
tical property is fully determined by the 4x4 correlation 
matrix x which is defined by 



- 1/ 



(3) 



where the vector x — [Xi , Yi , X2, Y2) and the quadrature 
operators Xj — {aj + aj)/\/2 and Yj = —i{aj — a])/-\/2. 
With the correlation matrix x, the continuous variable 
entanglement between the two cavity fields can be well 
quantified with the measure of logarithmic negativity. By 
re-expressing the correlation matrix x in terms of three 
2x2 matrices Qi, Q2, and ^3, 



X 



Qi 03 



the logarithmic negativity En{x) is defined as [3f 
i^jv(x) =inax[0,-ln(2A)], 

where 



A = 



A - - 4 Detx 
71 ' 



(4) 



(5) 



(6) 



(1) 



and A = Det^i -I- Det^2 — 2 Det^a. Thus, the entangle- 
ment between the two cavity fields occurs for i?Ar(x) > 0, 
i.e., for A > i. From the definitions of Eqs. ^ and (U), 
we see that the matrix elements of x ^-re a linear function 
of the second-order quantities 

n^Jit) = {alit)aJit)), = (a,(i)a,(t)) (7) 
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plus its hermitian conjugate. The entanglement dynam- 
ics between the two cavity fields at any time is then com- 
pletely determined by these time-dependent second-order 
quantities. 

On the other hand, the Hamiltonian of Eq. ([TJ de- 
scribes effectively the two spatially-separated nanocavi- 
ties coupled to a common reservoir with a controllable 
spectral structure. In other words, the waveguide plays 
a role of a structured reservoir as a mediator between 
two nanocavities. We can use the exact master equation 
we developed recently for nanocavities coupled to waveg- 
uides in photonic crystals [3TI-[33| to investigate the exact 
entanglement dynamics of the two nanocavities coupled 
to the waveguide in photonic crystals. By assuming that 
two nanocavities are initially uncorrelated to the struc- 
tured reservoir (i.e., the waveguide) and the waveguide 
is initially in vacuum, the exact master equation for the 
density operator p{t) of the cavity system can be ob- 
tained through the Feynman- Vernon influence functional 
approach 37[ in the framework of coherent state path- 
integral representation The resulting exact master 
equation is given by 



dpjt) 
dt 



-i[H,s{t),p{t)] 




(8) 



FIG. 2: (a) The spectral density Jaiio) for the cavities at 
different locations rii. (b) The spectral density Ji2(ci.') for 
the different sites n2 of the second cavity and the first cavity 
location ni = 1. Here the two cavity frequencies are taken 
to be equal and also equal to the resonator frequency in the 
waveguides: uj\ = uJ2 = uJo- The cavity- waveguide coupling 
^1 = ^2 = 0.2^0, and the hopping rate in the waveguide = 
0.05^0. 



Since the waveguide in photonic crystals has a narrow 
but continuous band structure, the spectral densities be- 
comes Jy(w) — 27r^?(cj)Vi(w)Vj (w), where p(a;) is the den- 
sity of states in the waveguide determined by Eq. 
Explicitly, we have 



where Hog{t) — ijJij(t)alaj is the effective Hamiltonian 
of the cavity fields with the time-dependent renormal- 
ized frequencies u}ij(t), which is resulted from the back- 
reaction of the the waveguide to the cavity system. The 
time-dependent coefficients ■j{t) describe the dissipation 
of the cavity system due to the coupling to the waveguide. 



The coefficients 
determined by 



UJ, 



j{t) and 7(t) are non-perturbatively 



(9a) 



l^J(t) = --m)^l-Ht)+h■C.],J, (9b) 

where fi{t) is the cavity photon propagating function 
which obeys the integrodifferential equation of motion 



d_ 
Jt' 



flit) = -lUflit) - f g{t- T)fl{T)dT, (10) 



with the initial condition fJ.(to) = 1. Here, the frequency 
matrix u> = Diag[wi,a;2] is a diagonal frequency matrix 
of the two cavities. The integral kernel in Eq. pH]) in- 
volves the time-correlation function gir — r'), which non- 
perturbatively characterizes the non-Markovian memory 
structure between the cavity system and the waveguide. 

By introducing the spectral densities of the waveguide: 
Jij{uj) = 2'K^j^VikVjk5{uj — ujk), the time-correlation 
function can be expressed as 



(11a) 



(12a) 



2e,sin[n,arcsin(^M^^H^)], (12b) 



where ujq — 2^o < < -I- 2^o is the band of the waveg- 
uide. In Fig. [21 the spectral densities Jii(w) — J22{^) 
and Ji2(i^) — J2i(s^) are plotted as a function of lo with 
different distances between the two cavities, for the case 
of the equal cavity frequencies loi = uj2 = ujc and the 
equal cavity- waveguide couplings = ^2 = ■Cc- As we 
will show later in the next section, the entanglement char- 
acteristics depend heavily on the spectral structures. In 
fact, it is the back-action through the off-diagonal ma- 
trix element Ji2(a;) of the spectral density that induces 
an effective coupling between the two separated nanocav- 
ities which leads to the entanglement generation. The 
environment-assisted continuous variable ent ang lement 
has been investigated in the literature [13, l40j |. How- 
ever, the entanglement in these investigations is uncon- 
trollable. Here, the waveguide (a structured reservoir)- 
induced entanglement between two spatially-separated 
nanocavities in photonic crystals is fully controllable and 
is promising for quantum information processing in all- 
optical circuits. 

Now, the exact temporal behavior of the second-order 
quantities of Eq. ([T]) can be completely determined by 
the exact master equation. Explicitly, from the exact 
master equation ([5]), it is not too difficult to find that 
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the second-order quantities obey the fonowing equations 
d 



dt 
d 

di 



n{t) = [fi{t)iJi-\t)]n{t) + n{t)[^l{t)^Ji-\t)]\ (13a) 



s{t) ^ [/i(t)M-^(t)]sW + s{t)[(i{t)i,-\t)Y . (13b) 



The exact sohitions of the above equations are found to 
be 

n{t) = t^it)n{0)fi\t) , sit) = /x(i)s(0)At^(t), (14) 

where n(0) and s{0) are the initial second-order quan- 
tities. Thus, once the photon propagating function fi(t) 
is solved from the integrodifferential equations of motion 
(ITU)) , the above exact solution allows us investigate the 
entanglement generation of the two separated nanocavi- 
ties. Meantime, the temporal evolution of the entangle- 
ment measure subjected to the non-Markovian dissipa- 
tion and fluctuation from the structured reservoir can be 
fully taken into account. 



III. RESULTS AND DISCUSSION 

Now we are ready to investigate the entanglement gen- 
eration between the two cavity fields and its temporal 
evolution for a given initially-separated state. We as- 
sume that the two cavity fields are initially prepared in 
a single-mode squeezed state, respectively, i.e.. 



|V'a.(0))=exp(-^a,2 



"farm 



i = l,2 (15) 



where the squeezing parameter ri is controllable as an 
input. The preparation of the initial squeezed states can 
be well accomplished by injecting into the nanocavities 
squeezed radiation fields produced via degenerate para- 
metric oscillation [4l[ . For the initial states of Eq. (ITSI) , 
the initial average photon numbers and two-photon cor- 
relations are given by 

nij{0) = sinh^(r,;)(5y , Sij(O) — sinh(ri) cosh(r,;)(5y . 

(16) 

With the help of the above initial conditions, we can 
analyze the temporal behavior of the correlation matrix 
X in Eq. and then the logarithmic negativity En 
given by Eq. ([5]). In the following, we will discuss the 
entanglement in three cases: (i) the cavity frequency uJc 
is resonant with the band center luq (wc = wq); (ii) the 
cavity frequency uJc stays the outside of the waveguide 
band (jwc — wo| > 2^o); and (iii) the cavity frequency lUc 
lies within the waveguide band (|wc — i^ol < 2^o)- 



At first, we consider the situation that the frequency 
LUc of the two cavities is resonant with the band center ujq 



of the waveguide. In this case, the steady-state average 
values in Eq. (jl4p in the weak cavity- waveguide coupling 
can be analytically obtained from Eq. (jA.3|) as 

t \ 722sinhV I 7iisinh2r 
(a[ai)s = ■ , (a^a2)s = — — — , (17a) 



7ll + 722 

/ 2\ _ 722 sinh2r ^ 2\ 



{aia2) 



2(711 + 722) ' 
7J2 sinh 2r 
2(711 + 722) '' 



711 + 722 
j'li sinh 2 r 

2(711 + 722) 



(17b) 



7i2 sinh r 
711 + 722 



(17c) 



where the subscript "s" denotes the steady state, 7^'^, = 
7m' 6"^*""=*, and the initial squeezing ri = r2 = r is 
also assumed. From the above results, the expression 
of the stationary logarithmic negativity E'^ can be ob- 
tained and the result is rather cumbersome so that we 
do not give it here explicitly. Nevertheless, it can be 
seen from Eq. p?| that the waveguide-induced collec- 
tive effect between the two cavities (characterized by the 
cross damping factor 712 in Eq. (jl7cp ) is crucial for gen- 
erating the entanglement. Since the cross damping rate 

712 = V711722 with 7i^ = (^)sin^(2^) in the weak 
coupling region, the stationary entanglement cannot be 
generated if any one (or both) of the sites ni and n2 
of the two cavities is even. Only when ni and n2 are 
both odd numbers, can the cross damping rate not van- 
ish. For the equal couplings = ^2 = ^c, we have 
7ii = 722 = 2^^/^o- Then, the steady-state logarithmic 
negativity Ef^ is reduced to 



E 



N 



(18) 



This shows that the entanglement degree between the 
two spatially-separated nanocavities can be controlled 
by the input initial squeezing parameter in the reso- 
nant case. Furthermore, from the definition of the pu- 



rity P = Tr(p^) = 



it is also not difficult 



4\/Dotx ' 

to find that the steady purity Pg = 1 for the steady 
state of Eq. ([TT]) . This can be clearly understood by 
performing an unbalanced beam-splitter transformation 
di = sin^fli + cos9a2 and ^2 — cos^oi — sm9a2 on the 
Hamiltonian in Eq. ([T]), where sin6' = Vik/\/V^^. + V^^.- 
Then the whole Hamiltonian of Eq. ([1]) is reduced to 
H = T.^'^cdU, + \/V^^+V^{dibl + bkd\). It shows 
that the waveguide only couples to the collective mode 
di , while the other collective mode ^2 is totally decoupled 
from the waveguide. Therefore, the collective mode d2 of 
the two cavities forms a decoherence-free subspace, which 
does not subject to any dissipation due to the presence of 
the structured reservoir. This results in the steady and 
pure entangled state between the two distant nanocavi- 
ties, as shown in Fig. |3] 

In Fig. [3l the exact temporal evolution of the loga- 
rithmic negativity Ej^ and the purity P of the two-mode 
cavity field are plotted. From it, we see that the station- 
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FIG. 3: The temporal behavior of the logarithmic negativity 
En and the entanglement state purity P for different values of 
the coupling r) and the site n2 the second cavity located. The 
cavity frequency ujc = i^o, the hopping rate = O.OStJo, the 
squeezing parameter r = 1.0, and the site of the first cavity 
ni = 1. 
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FIG. 4: The temporal behavior of the logarithmic negativity 
En (red lines) and the entanglement state purity P (blue 
lines) for difi'erent values of the coupling rj and the site n2 the 
second cavity located. The other parameters are the same as 
in Fig. 13 



~ 1 for all the cases in which the first cavity siting at 
ni = 1 and the second cavity siting at n2 = 5, 9, 15 with 
the coupling 77 = fc/Co = 0.08,0.05,0.01, respectively. 
It also shows that with the increase of the distance be- 
tween the two cavities, properly decreasing the cavity- 
waveguide coupling rj leads to the same steady and pure 
entangled state. This is because the oscillating profile 
of the spectral density Jii{uj), as shown in Fig. [21 be- 
comes narrower around the band center ujq as the site n2 
increases. This in turn requires longer time for achiev- 
ing the steady states as the coupling decreases. Thus, a 
proper choose of the cavity- waveguide coupling strength 
can yield the same entanglement degree and purity in the 
resonant case for the different distances between the two 
cavities, as shown in Fig. [H 

In Fig. |4l the entanglement and purity are plotted 
for the enhanced cavity- waveguide couplings r/ = 0.2 
and 0.4, with the different locations of the second cavity 
n2 = 5, 15 and 25. It shows that as the cavity- waveguide 
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FIG. 5: The temporal behavior of the logarithmic negativity 
En- The other parameters are the same as in Fig. [S] 



coupling increases, the entanglement decreases by accom- 
panying with some oscillations in the time evolution. The 
oscillation comes from the non-Markovian effect due to 
the back-action between the cavity and the waveguide 
when the cavity- waveguide coupling increases. Besides, 
one can also see that the entanglement sudden death and 
sudden birth occurs in the short-time regime with a rel- 
atively large cavity-waveguide coupling (77 = 0.4) and 
also a relatively long distance between the two cavities, 
see Fig. |H[b) and (c). In addition, as we see, the re- 
sulting entangled states become usually a mixed state 
as the cavity- waveguide coupling increases. By compar- 
ing Fig.|4](a)-(c), one can also find that the entanglement 
decreases as the increase of the distance between the cav- 
ities as well as the cavity-waveguide coupling. Therefore, 
to maintain the high entanglement for the two distant 
cavities, a small cavity- waveguide coupling is more fa- 
vorable in the resonant case. 

Furthermore, as shown in Fig. [5] the entanglement ap- 
pears in the case of the first cavity siting at rii = 1 and 
the second cavity siting at a even number of n2. The 
exact numerical result shows that the entanglement be- 
tween the two cavities can be generated under the rel- 
atively large coupling in a short time scale, but as the 
time goes the entanglement soon decays to zero. This is 
consistent with the analytical solution given by Eq. (|17l) , 
where it is pointed out that the steady-state entangle- 
ment in this case can not exist if any one (or both) of the 
sites ni and 712 of the two cavities is even. 



B. 



\luc - iOo] > 2^0 



To further investigate the controllability of the en- 
tanglement generation of the two spatially-separated 
nanocavities, we consider next the entanglement behavior 
for the cavity frequency lOc outside the waveguide band, 
i.e., jwc— wqI > 2^0- Let the first cavity site at ni = 1 and 
the second cavity locates at different sites. The temporal 
behaviors of the entanglement and the purity with differ- 
ent values of the cavity-waveguide coupling are plotted 
in Fig. [51 It shows that the entanglement exhibits a very 
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FIG. 6: The temporal behavior of the logarithmic negativity 
En and the purity P for the different values of the coupling 
rj and the site 722, while the cavity frequency ljc ~ 1.2ctJo 
which stays the outside of the waveguide band. The other 
parameters are the same as in Fig. [S] 



regular oscillation with a even stronger entanglement de- 
gree E'^'^^ > 1 for the input squeezing r — 1. For ex- 
ample, the maximal entanglement Ef^^^ « 1.94 for the 
coupling 7] = 0.2 in the present case, see Fig.lH^a). When 
the coupling is increased, the entanglement degree oscil- 
lates faster but the maximal entanglement is degraded a 
little bit only. 

The above phenomenon can be understood as follows: 
in the weak coupling region, the damping rates 7^^ = 
since the spectral densities Jij{LUc) = when the fre- 
quency uJc of the cavities lies outside the band of the 
waveguide. However, the cavity frequency shift, Suj = 
V f du-^-^^ does not vanish in this case. As a result, the 
reduced density matrix is purely determined by the effec- 
tive Hamiltonian i?eff = '^i{yJc + 5^jJii)o\ai + 5uji2{a\a2 + 
a^tti). In other words, an effective beam-splitter- type 
coupling (determined by the cross frequency shift 81012) 
between the two cavities is induced by the waveguide, 
which results in the entanglement for the separated 
squeezing inputs. With the weak coupling solution of the 
propagating function fibmit) given in the Appendix, we 
have explicitly yUbm.ii = cos(5wi2i)e~'"'^* and = 
— i sin(5cLii2te~*'^'=*, where Cjc = i^c + and SUJ22 ~ 
Sujii. The average values in Eq. reduce to (a|ai) = 
sinh^ r, (af) — sinhr coshr cos(2(5wi2t)e~^''^'^*, (0102) = 
— I sinhr coshr sin(25a;i2i)e~^*"'^*, and {a\a2) = 0. Thus, 
unlike the resonant case, the entanglement in this case is 
purely determined by the cross frequency shift Suji2. At 
the times when cos(2(5a;i2<) = 0, we have the nonzero av- 
erage values (a|ai) ~ sinh^ r and |(aia2)| — sinhr coshr, 
which corresponds to a pure two-mode squeezed vacuum 
state with the squeezing parameter r |42j. Accordingly, 
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FIG. 7: The temporal behavior of the logarithmic negativ- 
ity En and the purity P for the cavity frequency within the 
waveguide band. Here we take Uc — l.OSctJo and the site 
n2 — 5. The insert in the figure depicts the average photon 
numbers {a\ai) of the two cavities. The other parameters are 
the same as in Fig. [H] 

the entanglement degree at these times becomes optimal 
with the maximal logarithmic negativity 

EY^""" = 2r. (19) 

Fig. |6ja) shows that the maximal entanglement ap- 
proaches the limit 2r when the coupling becomes suf- 
ficiently weak and the weak coupling solution given in 
Appendix is almost exact. With the cavity- waveguide 
coupling increase, the maximal entanglement is decreased 
a little bit, as shown in Fig. IH^a). 

Fig.lHKb) depicts the corresponding purity of the entan- 
gled state. Since the damping rate vanishes in the weak 
coupling limit when the frequency of the cavities lies out- 
side the waveguide band, the entangled state should be 
a pure state. Fig. HKb) shows that for rj — 0.2, the ex- 
act numerical solution gives the purity P « 0.97 which 
is consistent with the weak coupling solution. When 
the cavity-waveguide coupling is increased, the purity 
of the waveguide-generated entanglement is decreased. 
Fig. [6Kc) and (d) show further the dependence of the 
entanglement degree and the purity on the distance be- 
tween the two nanocavities. Similar behavior of the os- 
cillating entanglement is obtained as the cavity distance 
changes. This is due to the cross frequency shift 5uji2 
again, which is determined by the cross spectral density 
Ji2(a;) that decreases slightly as n2 increases for the fixed 
ni = 1. Besides, it is also found that the maximal entan- 
glement and purity in the long-time region do not change 
obviously with the changes of the distance between the 
two nanocavities. 



C. \uJc -L0o \ < 2^0 

Finally, we shall consider the case that the cavity fre- 
quency UJc is not resonant with the band center but it 
still lies inside the waveguide band. In Fig. 7 (a), the en- 
tanglement and the purity are plotted for the cavity sites 
ni = 1 and n2 = 5 with the cavity frequency uJc — l.OSwo 
at which the spectral density J22(w) = for uj = uJc 
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FIG. 8: The temporal behavior of the logarithmic negativity 
En for the cavity frequency ujc = 1.06a;o and the second cav- 
ity locating at the site n2 = 5. The other parameters are the 
same as in Fig. O 

(see the red-dashed line in Fig. [2]). In this case, the phe- 
nomenon of entanglement sudden death and sudden birth 
occurs in the weak coupling region. The existence of the 
entanglement sudden death and sudden birth originates 
from the fact that the damping 722 ~ and also 712 ~ 
for ojc — l.OSwo so that the entanglement is purely gov- 
erned by the cross frequency shift 5uJi2, which leads to 
the entanglement oscillation. In the long time limit, the 
entanglement is destroyed through the damping channel 
711 0. When the cavity-waveguide coupling increases, 
the entanglement only exists in a very short time and 
then quickly decay to zero significantly, which is quite dif- 
ferent from the resonant case shown in Fig. 21 The insert 
in Fig. El^a) is the intracavity average photon numbers 
{ctjCij) which approach to zero in the long-time limit. It 
tells that the cavity evolves asymptotically into a vacuum 
state. This is why the entanglement in the steady-state 
vanishes. The purity shown in Fig. El^b) approaches to 
one in the steady-state limit because the steady state is 
just a trivial vacuum state. 

In Fig. [HI we plot the entanglement for the cavity fre- 
quency LUc = l.OGwo, while the spectral density Jaioj) ^ 
(z = 1, 2) at cj = LUc- Then, both the decay channels 
711 and 722 become active, and the entanglement is a 
combination effect of the the non-zero cross damping 712 
and the non-zero cross frequency shift 5u)i2. Compared 
to that in Fig. [TJa), the entanglement here is enhanced 
significantly in the long-time regime. The phenomenon 
of the entanglement sudden death and sudden birth dis- 
appears in this case. As a result, we see that the entan- 
glement generation can be easily controlled by changing 
the cavity frequency within the band of the waveguide. 



IV. CONCLUSIONS 

In conclusions, the generation of continuous variable 
entanglement between two spatially-separated nanocavi- 
ties mediated by a CROW in photonic crystals is inves- 



tigated. By solving the exact master equation for the 
cavity system coupled to the waveguide as a structured 
reservoir, the entanglement and the purity of the two- 
mode cavity field is analyzed in details for the initially- 
separated squeezing inputs. We found that the steady 
and pure entangled state of the two distant cavities can 
be generated for a weak cavity-waveguide coupling when 
the cavities is resonant with the band center of the waveg- 
uide. It also shows that the strong coupling of the cav- 
ities to the waveguide can lead to the phenomenon of 
entanglement sudden death and sudden birth. When the 
cavity frequencies lie outside the waveguide band, the 
optimal entanglement can be achieved by the cross fre- 
quency shift between the two cavities, which is induced 
by the waveguide. When the cavity frequencies are not 
resonant with the band center but still within the waveg- 
uide band, the entanglement can also exhibit the phe- 
nomenon of sudden death and sudden birth even for a 
weak cavity-waveguide coupling. By changing the cavity 
frequencies within the band of the waveguide, the en- 
tanglement and the occurrence of entanglement sudden 
death and birth can be easily controlled. These inter- 
esting results show that the waveguide, as a controllable 
structured reservoir, can be used to entangle efficiently 
the two spatially-separated nanocavities. More impor- 
tantly, only a weak cavity- waveguide coupling is required 
for achieving the optimal entanglement. We expect that 
these interesting features can be realized in experiments 
and find its application in continuous variable quantum 
communication networks. 
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Appendix: A weak-coupling analytical solution 

It should be pointed out that the master equation of 
Eq. ([H]) is exact, far beyond the Born-Markovian approxi- 
mation and valid for arbitrary cavity- waveguide coupling. 
The back-action between the cavities and the waveguide 
is embedded into the time-dependent coefficients uJij(t) 
and 7(t) of Eq. which are in turn determined com- 
pletely by the propagating function fi{t) of Eq. (flO)) . 
Generally, it is not easy to obtain the analytical prop- 
agating function /x(t). However, for a weak coupling, the 
analytical solution of the photon propagating function 
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can be found as 



weak-coupling limit is then given by 



(A.l) 



where the damping rate becomes time-independent: 7 = 
J{lUc)/2, and the waveguide-induced renormalized fre- 
quency is given by a; = Ucl + Sui with 



Suj =V 



duj J{lo) 

l-K UJ ~ Ijj 



(A.2) 



Here V denotes the principle value of the integral. When 
the cavity frequency uJc is resonant with the band center 
of the waveguide (wc = wq), the frequency shift — 
0, and the explicit propagating function /Xf„„(t) in the 



1 



Mbm.ll 



711 + 722 



[722 +7iie"^^"+^^^^*]e-'"=*, (A.3a) 



Mbm,22 = ^ [711 + 722e-(^"+^-)*]e-''^=*, (A.3b) 

711 + 722 



Mbm,12 — A* bm,21 — 



712 



711 + 722 



_ p-(7ii+722)il„-icJct 



[1-e 



(A.3c) 



where the damping rates 7ji = (^)sin^(^^) and the 
collective damping 712 = ^711722- 
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